SUBSIGNATURES OF SYSTEMS 



JEAN-LUC MARICHAL 



Abstract. We introduce the concept of subsignature for semicoherent sys- 
tems as a class of indexes that range from the system signature to the Barlow- 
Proschan importance index. Specifically, given a nonempty subset M of the 
set of components of a system, we define the M-signature of the system as the 
|M|-tuple whose fcth coordinate is the probability that the fcth failure among 
the components in M causes the system to fail. We give various explicit linear 
expressions for this probability in terms of the structure function and the dis- 
tribution of the component lifetimes. We also examine the case of exchangeable 
lifetimes and the special case when the lifetime are i.i.d. and M is a modular 
set. 



1. Introduction 

Let S = (C,cj),F) be an n-component system, where C = {1, . . . , n} denotes the 
set of components, qb denotes the associated structure function </>:{0,l}" -»■ {0,1} 
(which expresses the state of the system in terms of the states of its components), 
and F denotes the joint c.d.f. of the component lifetimes T\,... ,T n . We assume 
that the system is semicoherent, i.e., the structure function <f> is nondecreasing in 
each variable and satisfies the conditions 0(0, . . . , 0) = and 0(1, . . . , 1) = 1. We 
also assume that the c.d.f. F has no ties, that is, Pr(Ti = Tj) = for all distinct 
i,j e C. 

The signature of the system, a concept introduced first in 1985 by Samaniego [9] 
for systems whose components have continuous and i.i.d. lifetimes and then recently 
extended to non-i.i.d. lifetimes (see [5] and the references therein), is defined as the 
n-tuple p = {pi, . . . ,p n ), where pk is the probability that the fcth component failure 
causes the system to fail. That is, 

Pk = Pr(T 5 = T fe: „), fce{l,...,n}, 

where Ts denotes the system lifetime and Tk.n denotes the fcth smallest lifetime, i.e., 
the fcth order statistic obtained by rearranging the variables T\, . . . , T n in ascending 
order of magnitude. 

In the i.i.d. case (and even when the component lifetimes are exchangeable) the 
system signature p is independent of the c.d.f. F . It is then called the structural 
signature and denoted by s = (si, . . . , s n ), where Sk = Pr(Ts = Tk-. n ). In fact, Boland 
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[3] showed that Sk can be explicitly written in the form 

(1) s k = Z Ax<K*)- E At HA). 

AzC \\ A \) AzC {\ A \) 

|A|=n-fc+l \A\=n-k 

Here and throughout we identify Boolean vectors x 6 {0,1}™ and subsets A Q C by 
setting Xi = 1 if and only if i € A. We thus use the same symbol to denote both a 
function /:{0,1}™ -*■ R and its corresponding set function f:2 c -*■ K, interchange- 
ably. For instance, we have (j>(C) = 0(1, . . . , 1). 

Marichal and Mathonet [5] extended formula (JTJ to the non-i.i.d. case into 

(2) Pk = E l(A)HA)- E Q(A)HA), 

A£C AsC 
|j4|=n-fe+l |A|=n-fc 

where the function q:2 -> [0, 1], called the relative quality function associated with 
F, is defined by 

q(A) = Pr I max Tj < minT^ ) . 

That is, q(A) is the probability that the best \A\ components are precisely those in 
A. 

Since its introduction the concept of signature proved to be a very useful tool 
in the analysis of semicoherent systems, especially for the comparison of different 
system designs and the computation of the system reliability (see, e.g., [TU] for the 
i.i.d. case and 0E] for the non-i.i.d. case). 

The Barlow-Proschan importance index of the system, another useful concept 
introduced first in 1975 by Barlow and Proschan pQ for systems whose components 
have continuous and independent lifetimes and then extended to the general de- 
pendent case in [H H], is defined as the n-tuple Ibp whose jth coordinate is the 
probability that the failure of component j causes the system to fail, that is, 

= Pr(T s = Tj) . 

Just as for the signature, in the i.i.d. case (and even when the component lifetimes 
are exchangeable) the index Ibp is also independent of the c.d.f. F. It is then called 
the structural importance index and denoted b = (&i, . . . ,b n ), where bj = Pr(Tg = 
Tj ) . An explicit expression for bj in terms of the structure function values is given 

by 

(3) ^ = £ _L-A^(A), 

AzCs{j} n {\A\ ) 

where Aj4>(A) = <p{Au {j}) - <j>(A \ {j}). Marichal and Mathonet [B] extended this 
formula to the non-i.i.d. case into 

(4) 4p = E <Z,-U)A^(A), 

Ac C s{j} 

where the function qj\2 c ^' -*■ [0,1] is defined by 

(5) qJA) = Pr max Ti = T,-< min Ti ] . 

J V i t C\A J if A I 

That is, qj(A) is the probability that the components that are better than compo- 
nent j are precisely those in A. 



3 



For instance, when n = 4 we have 

<&({1,3}) = Pr(T 4 <r 2 <min{T 1 ,T 3 }) 

= Pr(T 4 < T 2 < Ti < T 3 ) + Pr(T 4 < T 2 < T 3 < T t ) . 
More generally, ([5]) can be rewritten as 

(6) q 3 (A) = Y, Pv(T aW <-<T <n) ), 

(T£e„ :{o{n-\A\+l),...,(j{n)}=A 
<r(n~\A\)=j 

where & n denotes the set of permutations on C. 

These important concepts of signature and Barlow-Proschan index motivate the 
introduction of the following more general concept. Let M be a nonempty subset 
of the set C of components and let m = \M\. We define the M -signature of the 
system as the m-tuple p M = (pf 1 , ■ ■ . ,p*'), where p^ 1 is the probability that the 
kth failure among the components in M causes the system to fail. That is, 

pf = Pr(T s = T k:M ), ke{l,...,m}, 

where T^-.m denotes the fcth smallest lifetime of the components in M, i.e., the kth 
order statistic obtained by rearranging the variables Ti (i e M) in ascending order 
of magnitude. A subsignature of the system is an M-signature for some M £ C. 
We naturally denote the probability p^ by s^f whenever the component lifetimes 
are exchangeable. 

Clearly, when M = C the M-signature reduces to the signature p, which shows 
that the signature is a particular subsignature. At the opposite, when M is a single- 
ton {j} the M-signature reduces to the 1-tuple = (p\^), where p\^ = Pr(Ts = 
Tj) is the jth coordinate of the Barlow-Proschan index. Thus, the subsignatures 
define a class of 2 n - 1 indexes that range from the standard signature (when M = C) 
to the Barlow-Proschan index (when M consists of a single component). 

Remark 1. The concept of M-signature is particularly relevant for instance when 
the system has a huge number of components and M is a small subset of potentially 
unreliable components. 

Interestingly, we have the following link between the subsignatures and the 
Barlow-Proschan index. For every nonempty subset M £ C, we have 

m 

Y, Pk = Pr(T s = T, for some j e M) = £ /<•» . 

fc=l jeM 

If this probability is not zero, then we can express the normalized M-signature 
Pk 1 1 SSi pT as t ne conditional probability 

Jn M = Pr(T s = T k:M | T s = Tj for some j e M) . 
Le=i Pi 

Example 1. Consider the six-component system indicated in Figure [TJ where the 
components have i.i.d. lifetimes. The corresponding structure function is given by 

(f)(xi,...,Xe) = (xiUX2)x 3 X4,(x 5 UX 6 ) , 

where u is the coproduct operation defined by x u y = 1- (1 — x)(l - y). 

Table [TJ gives the normalized subsignatures for some subsets (those for the re- 
maining subsets can be derived by symmetry). These normalized subsignatures 
illustrate the following fact: The more concentrated the distribution of s M is on 
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M 




{1,2} 


(0.000, 


1.000) 


{1,3} 


(0.846, 


0.154) 


{3,4} 


(1.000, 


0.000) 


{1,2,3} 


(0.467, 


0.533, 0.000) 


{1,3,4} 


(0.833, 


0.167, 0.000) 


{1,2,5} 


(0.167, 


0.666, 0.167) 


{1,2,3,4} 


(0.538, 


0.462, 0.000, 0.000) 


{1,2,5,6} 


(0.000, 


0.500, 0.500, 0.000) 


{1,2,4,5} 


(0.412, 


0.471, 0.118, 0.000) 


{1,3,4,5} 


(0.692, 


0.269, 0.038, 0.000) 


{1,2,3,4,5} 


(0.464, 


0.429, 0.107, 0.000, 0.000) 


{1,2,3,5,6} 


(0.263, 


0.421, 0.316, 0.000, 0.000) 


{1,2,3,4,5,6} 


(0.333, 


0.400, 0.267, 0.000, 0.000, 0.000) 



Table 1 . Normalized subsignatures for representative subsets 



the left, the more M is "critical" in the sense that each of its components makes 
the system vulnerable and potentially unreliable. □ 
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FIGURE 1. A six-component system 

In this paper we give explicit linear expressions for M-signatures, including a 
simultaneous generalization of formulas ([2]) and (|4]) as well as a linear expression 
in terms of the Mobius transform of the structure function (Section 2). We also 
investigate the special case when the component lifetimes are exchangeable (Section 
3). Finally, we examine the case when the lifetimes are i.i.d. and M is a modular 
set (Section 4). 

2. Explicit expressions 

In this section we provide various explicit expressions for the probability p^f = 
Pr(Ts = T/c;m). We start with a linear expression in terms of the functions qj and 
Aj</>, thus generalizing formula 

Theorem 2. For every nonempty set M £ C and every k e {1,. .. ,m}, we have 
(7) Pf = E E A 3 4>(A) = E E QjW^tiA). 

A^C jeM\A juM AsC\{j} 

\M\A\=k \M\A\=k 

Proof. Fix a e & n , set j £ C such that Tj = Tk-.M, an d let j' = a~ l {j). Then we have 
Pr(T s = Tj | T CT(1) < .« < T ff(n) ) = 0({a(j'), ■ • ■ " 4>{W + !),•■ • M")}) ■ 

Indeed, the left-hand expression of this equation takes its values in {0, 1} and is 
exactly 1 if and only if {a(l), . . . , 1)} is not a cut set and {c(l), ■ ■ ■ , cr(j')} is 
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a cut set. (Recall that a subset K £ C of components is a cut set for the function 
(f> if 0(C \ if) = 0.) That is, 

cf>({a( f ),..., *(n)}) = l and </>{{*(/+ l),...,<r(n)}) = 0. 

Since j' depends on tr, we denote it by j^.. By using the law of total probability, 
we then obtain 

Pr(T s = TW) = £ (^(Mi;),...,a(n)})-^(Wj;+l),...,<7(n)})) Pr(T ff(1) < - <T ff(n) ). 

Grouping the terms for which {<j{j' a ) 1 . . . , <r(n)} is a union Bi u B 2 u {i} of fixed 
sets, where i = cr(j^), Bi £ C\M, and £? 2 £ M\{i}, \B 2 \ = m-k, and then summing 
over these sets, we obtain 

Pr(T 5 = T fc:M ) E E E (0(S 1 uS 2 u{ i })-0(S 1 uS 2 )) 

1-^2 \=m-k 

E Pr(T CT(1) <-<r a(n) ). 

<t<eS„ :{tr(n-|S 1 uB 3 [+l),--,o-(n)}=BiUS 2 
o-(ri-|BiuB2|)=i 

By ([6]) the inner sum is precisely qi(Bi u B 2 ). We then conclude by setting A = 
£?i u B 2 . The second expression can be obtained by permuting the sums of the first 
expression. □ 

Example 3. Consider a 3-component system whose structure function is given by 

4>(xi,x 2 ,x 3 ) = (xiux 2 )x 3 = XlX 3 + x 2 x 3 - X\X 2 X 3 . 
For such a system, we have for instance 
P?' 3} = E E V(A)A^(A) 

As{l,2,3} je{l,3}\A 
|{1,3}sA|=1 

= ?i({3}) + (fe({l}) + fe({l,2}) 

= Pr(T 2 < Ti < T 3 ) + Pr(T 2 <T 3 < 1\) + Pr(T 3 < T x < T 2 ) + Pr(T 3 < T 2 < Ti). 
Interestingly, Theorem [5] provides new expressions for p^, namely 

Pk = E E Aj 4>(A) = E E QM)^<t>(A). 

AzC j"eC\A jeC A^C\{j} 

\A\=n-k \A\=n-k 

Replacing Aj <fi(A) by cf)(A u {j}) - <fi(A) in formula (JT)), we derive the following 
alternative expression for the probability p^f . 

Corollary 4. For every nonempty set M £ C and every fee {1, . . . , m}, we have 
(8) Pf = E E (-IJ^^AxW)^). 

jeM AcC 

|(MxA)u{i}|=/= 

Proof. The right-hand side of (O can be written as 

E E ij(A^{j})<t>(A)- e E <&U)^) 

jeM AcC, Aaj jeM AsC, Aj!j 

|A/\A|=fc-l \M\A\=k 

= e E ?,-(A)^u{i})- E E 

jeM AsC\{j} 3'eM A£C\{j} 

|M\A|=fe |A/\A|=fc 
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which is precisely the right-hand side of ((TJ). □ 

We now provide an alternative linear expression for the probability pjf that 
generalizes formula ([2]). For every nonempty set M £ C, define the set functions 
q l M : 2 C \ {0} -> K and qj^ : 2 C \ {C} -+Rby 

9m ( A ) = E 9i ( ^4 n {i} ) = Pr (3 j 6 M : max T, < 7) = min ) 



and 



9m (^) = E ^ M) = Pr (3 i 6 M : max = T s < min I) ) 

jeM\A «OvA itA 



respectively. 

Corollary 5. For every nonempty set M £ C and every k e {1, . . . , m}, we have 

(9) pff = E q {M)^A)- y ilM)^(A). 

AsC AeC 
\MnA\=m-k+l \MnA\=m-k 

Proof. By ([2|) we have 

Ptf = E E V(A)A^(A) 

AsC jzM\A 
\M\A\=k 

= E E qj(A)<f>(Au{j})- £ E lAA)4>(A) 

AsC jtM\A AsC jtM\A 

\M\A\=k \M\A\=k 

E E V(A-{J})HA)- Y Qm(A)HA), 

A^C jeMnA A^C 
|M\A|=fc-l |M\A|=fc 

which completes the proof. □ 



We end this section by providing an explicit linear expression for the probability 
p^f in terms of the Mobius transform of the structure function <f>. Recall that the 
Mobius transform of <f> is the set function c^: 2 C -*■ K which gives the coefficients of 
the unique multilinear expression of the structure function, that is, 

0w = e c M) n^- 

AeC ieA 

The conversion formulas between and <fr are given by 

(10) C</> (A) = ^(-1)I A HB| 0(S) and HA) = Yc^B). 

B^A B^A 

A very simple linear expression for p^f in terms of the Mobius transform is given 
in the following theorem. 

Theorem 6. For every nonempty set M £ C and every k e {1, . . . ,m}, we have 

pi 1 = Y C 4>(A) Pr(T fe:A f = minTi) , 

aTc V teA ' 

\MnA\im-k+l 

or equivalently, 

pf = Y c^A)Pr(T kM = mmT l ). 
aTc v ieA ' 
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Proof. By substituting the second formula of (fT0|) in (|9|) and then permuting the 
resulting sums, we obtain 

Pk = E q' M {A) Y,m{B)~ £ q ' M (A) £ m(B) 

AoC B^A AoC B^A 

\MnA\=m-k+l \MnA\=s-k 

E E <4U)- E E qK a )- 

B^C A?B B^C A^B 

\MnB\iin-k+l AfnA|=m-fc+l \MnB\im-k \MnA\=m-k 

However, we have 

Y, 1m{A) = PrhjtM, 3 A 2 B, \MnA\ =m-k+ 1: max Ti<Tj = minT;) 

A^B 
\MnA\=m-k+l 



= Pr [T k :M ^ minT; ) 
E Qm( A ) = Pr ( 3 j e M, 3A^B, \M n A\ = m - k : max T = Tj < mm T j 
= Pr(T fe:M <minT i ). 



and 



A2B 
\MnA\=m-k 



Thus, we have 

Pk = E m( J B)Pr(T fe:M <minTi)- £ m(S) Pr(T fe:M < minT 4 ) 

\MnB\im-k+l \MnB\im-k 

Y, m(B)Pr(T k .. M $mmT z )+ £ m(B) Pr (T fc:M = minT,) . 

\MnB\=m-k+l \MnB\im-k 

We observe that we cannot have T^.m < min^s T if |MnS| = m-k+1. This proves 
the first formula of the theorem. To see that the second formula holds, just observe 
that we cannot have Tk-.M = niin^s T if \M n B\>m—k+l. □ 

Example 7. Consider a system defined as in Example 02 that is, with a design 
given by the (multilinear) structure function 

4>(xi,X2,X 3 ) = X1X3 + x 2 x 3 - XiX 2 X 3 . 

For such a system, we have for instance 



p\ 1,3} = E ^(A)Pr(T 1:{1)3} =mhiT i ) 

As{l,2,3} KA 



1 + Pr ( min T = min T J - Pr ( min T = min T ) 

W{1,3} ie{2,3} ' Me{l,3} ie{l,2,3} ' 

= Pr(T 2 < T x < T 3 ) + Pr(T 2 < T 3 < r x ) + Pr(T 3 < T x < T 2 ) + Pr(T 3 < T 2 < Ti). 

Remark 2. We observe that the probability Pr(Tk-,M = rnin^ T) is exactly the fcth 
coordinate of the M-signature of the semicoherent system obtained from the current 
system by transforming the structure function into <?i(x) = YliaA x i- This result 
follows immediately from the fact that the modified system has lifetime min^Ti. 

From Theorem [5J we immediately derive the following corollary, which was al- 
ready established in [5J. 
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Corollary 8. For every k € {1, . . . , n} and every j e C , we have 
Pk = E c <t>(A) Pr(Tfe :n = minTA and J^p = E c </>(^) Pr (^j = minTi) . 

3. Exchangeable component lifetimes 

We now consider the special case when the functions qj (j e C) satisfy the 
condition 

(11) Qj(A) = -t^, Ac C s{j}, jeC. 

It is easy to see that this condition holds whenever the lifetimes Tj., . . . ,T n are i.i.d. 
or, more generally, exchangeable (see |6]). In this case, for every nonempty subset 
M £ C, we also have 



.A/ 



and 



Theorem [2] and Corollaries 2] and [5] are then immediately specialized to the 
following result. Here the probability pf? is naturally denoted by s^f . 

Corollary 9. Assume that the functions qj (j e C) satisfy condition ill]) . For 
every nonempty set M £ C and every k € {1, ... ,m}, we have 

(12) *f = E ""T^K E A,- 0(A) = X E A, 0(A) , 

|M\A|=fc [Af\A|=fc 

E E (~^ A -A^HA) , 

|(MxA)u{j}|=fc 

4 1 - E ! ^^)- E -4^ HA). 

A^C n \\A\-l) AzC n \\A\) 

\MnA\=m-k+X \MnA\=m-k 

From (|12p we immediately derive new expressions for Sk, namely 

* = E -4^ E ^HA) = E E -tLa^a). 

| A | J jeCxA jeC AcCx{j} "A |A| J 

|A|=n-fc |A|=n-fc 

An expression for s^f in terms of the Mobius transform of the structure function 
is given in the following result. 

Corollary 10. Assume that the functions qj (j e C) satisfy condition hll\). For 
every nonempty set M £ C and every k e {1, ... , m}, we have 

m-\M,A\ 



S k ~ E C <t>( A ) , , . Ws + 



J. l\M\A\ + \A\\ 

A^C K { k ) 

k-li\M\A\<im-l 



or equivalently, 



fc /|MxA|+|Ah 



r fc' ' ) 
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Proof. For every A £ C, let </>a(x) = YlieA x i- F° r every B £ C, we then have 
Aj </>a(-B) = 1, if j e A and A s {j} E B, and 0, otherwise. 
Combining Remark [5] with (|T^|) , we then obtain 

Pr(r fc:M = mmr,) = £ £ _La^ a (B) 

jMsB|=fc 

E E 



jcMnA ix{j}c B cC\{ 3 } n C|B|) 
|M\B|=fc 

Partitioning ^4 into A\ = An M and A2 = A \ M, the latter expression becomes 
EE E 7 77^1 \ 

I -Si |=ti— k 

- ,L[ m _ k _ lAil + 1 ) b L M { h _ lA2l )J li -*J 

1*1 W) 



n+ iA 2 | (^i- 1 ) ■ 

We then conclude by Theorem [5] □ 
From Corollary 1101 we immediately derive the following result. 

Corollary 11. Assume that the functions qj (j e C) satisfy condition For 
every k e {1, . . . ,n} and every j e C , we have 

\A\ M Al ) 1 
s k = £ c,{A)^^-L and b . = £ ^(A) — . 

AzG K \ k ) AsC,ABj I 71 ! 

|j4|<n-fc+l 



4. SUBSIGNATURES ASSOCIATED WITH MODULAR SETS 

Suppose that the component lifetimes are i.i.d. and that the system contains a 
module (M, x), where M is the corresponding modular set and x- {0> 1} M "* {0, 1} 
is the corresponding structure function. In this case the structure function of the 
system expresses through the composition 

(13) 0(x) = ^( X (x M ),x^ M ), 

where x M = (xJum, x CxM = (xi) ieC ^M, and, setting C M = (C \ M) u {[M]}, 
the function ip: {0, 1} Cm {0, 1} is the structure function of the reduced system (of 
n—m+l components) obtained from the original system by considering the modular 
set M as a single macro-component [M]. For general background on modules, see 
Chap. 1]. 

Denote by Tm the lifetime of the module. By (|12[) the signature of the module 
is an m-tuple whose fcth coordinate is given by 

(14) Pr(T M = T fcM ) = £ Y. A 3 X(A). 

ASM m \\A\) i*M\A 
\A\=m-k 
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Denote by h^: [0, \\ Cm -*■ R the reliability function of the structure function t/j, that 
is, the multilinear polynomial function 

Mx) = e ^) n 

A^Cm it A ieCM^A 

We can easily see that the average value over [0, 1] of the function 

(15) (d [M ]M (*,•••,*) = E * |A| (i-*r- m - |A| A [M] ^) 

A£C M s{[Af]} 

(see [5] for details) with respect to the uniform distribution is given by 

[\d [M] h^)(t,...,t)dt = £ 1 ^ A [M] ^(A). 

Jo ac Cm n{[a/]} (n-m + l)( |A| J 

By (|3]) this quantity is precisely the [M]-coordinate of the Barlow-Proschan impor- 
tance index for the reduced system (Cm,^)- Thus, we have 

(16) j\d [M] h 4 ,)(t,...,t)dt = Pv(T s = T M ). 

The following result shows that s¥ factorizes into the product of Pr(Tj\/ = T^.m) 
and the expected value of the function {d[M]hip){t 1 . . . ,t) with respect to a beta 
distribution. When M is a singleton, this result reduces to (|16p . 

Theorem 12. Assume that the components have i.i.d. lifetimes. For every nonempty 
modular set M £ C and every fee {1, . . . , m}, we have 

r 1 

si' = Pr(T M = T fc:M ) J o r k ,m(t)(d[ M ]hii 1 )(t,...,t)dt, 

where ip is the structure function defined in il.3\) and rk, m (t) is the p.d.f. of the 
beta distribution on [0, 1] with parameters a = m - k + 1 and f3 = k. 

Proof. Let [0, l] c -»■ K and h x : [0, 1] M -+ R be the reliability functions of 4> and 
X, respectively. We then have 

Mx) = ^(/i x (x M ),x CsM ). 

By the chain rule it follows that, for every j € M, 

c^(x) = (5 [M] ^)(/ lx (x M ),x^ M )5 J / lx (x M ). 

Since d[M]h^ does not depend on its [M]-variable, we have 

(s [M ]M(Mx M ),x CxM ) = (d [M] h^)(0,^ M ). 

Since the reliability functions are multilinear, for every A\ £ M and every Ai £ 
C \ M we immediately obtain the factorization 

A j </>(A 1 uA 2 ) = A [M] iP(A 2 )A jX (Ai). 
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Therefore, by ([12]) and {H} we obtain 



S k EE/ n-l \ 

|Ai|=m-fc 



T E A 



,j<l>(Ai uA 2 ) 



A 2 £C\M 

and we conclude by (|T5|) . 



= Pr(T M = T fcM ) E 



|Ai|=m- 



E 



k 



m 




A [M] V(A 2 ) 



□ 



Interestingly, Theorem [T2] shows that in the i.i.d. case the ratio s¥ / Pr (Tm = 
Tu-.m) is independent of the structure x of the module (assuming Pr(Tj^ = T^.m) + 
0). Thus, if the components of the module are reorganized so that \ is modified 
into x' an d if the probability Pi(Tm = T^-m) remains unchanged, then so does s^f . 
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